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Abstract: Lattice artefacts are used, through modified lattice actions, as a 
tool to find the largest instantons in a toroidal geometry [0,L] 3 x [0, T] for 

£T) '. T — > oo. It is conjectured that the largest instanton is associated with 

tunnelling through a sphaleron. Existence of instantons with at least 8 pa- 
rameters can be proven with the help of twisted boundary conditions in the 

'"" I" 1 ■ time direction. Numerical results for SU(2) gauge theory obtained by cooling 

^h. are presented to demonstrate the viability of the method. 

: 

'. 1 Introduction 

X' 

Since the time of the discovery of instantons ffl in non-abelian gauge theories, as vacuum 
to vacuum quantum mechanical tunnelling events , their role in strongly interacting the- 
ories has been controversial, both in the continuum |§ and in the lattice formulation Q. 
For the continuum this has been mainly due to applying semiclassical techniques, which 
cannot be justified at strong coupling. In the lattice formulation the main problems were 
the instantons localised at the scale of the lattice cut-off for which topological charge can- 
not be defined unambiguously ||, and which have actions considerably lower than the 
continuum action of 8n 2 . Strictly speaking, there are no locally stable solutions on a lat- 
tice using the standard Wilson action ||, because this lattice action decreases when the 
instanton becomes more localized [[/J, as we will demonstrate also from analytic consid- 
erations. On a trial and error basis, different (improved) lattice actions were considered, 
some of them indeed giving rise to stable lattice solutions |§ . This paper will provide the 
proper framework to understand the stability. 

It is not too difficult to understand the reason of the instability. At finite lattice sizes the 
lattice action deviates from the continuum and this deviation is larger for stronger fields. 
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For the Wilson action, as we will show, the lattice artefacts make the action decrease as 
compared to the continuum. In the continuum, instantons have a scale (or size) parameter 
p, on which the action does not depend. But the smaller p becomes, the larger the fields 
get, which makes the lattice action decrease . On dimensional grounds one easily argues 
that (generically) Si att (a, p) = 87r 2 (l + (a/p) 2 d 2 + 0(a/p) 4 ) for p ^> a, which will be 
demonstrated in more detail further on. For the Wilson action H d% < 0, explaining the 
instability. Hence, one simply modifies the action, such that d 2 > 0, in order to get stable 
solutions for the maximal value of p allowed by the volume [0, L] 3 , which is kept finite. As 
we are interested in the classical solutions to the equations of motion, the modified action 
need not be of the type of an improved action || , for which typically one wants to achieve 
d 2 = 0, as in that case (as we will show) the (a/p) 4 term might still destabilize the solution. 

We deliberately want to keep d 2 > 0, which we will hence call over-improvement. The 
reason is, that our motivation for embarking on this project was to find the instantons 
with the largest scale p. This presumably will correspond to tunnelling over the lowest 
energy barrier, separating two classical vacua. The configuration that corresponds to the 
lowest barrier height is then conjectured to be a sphaleron (which exists due to the fact 
that we keep the volume finite). A sphaleron |TD[ is by definition a saddle point of the 
energy functional with precisely one unstable direction, which corresponds to the direction 
of tunnelling. In this way we use the instantons to map out the part of the energy functional 
relevant for the dynamical region where a semiclassical analysis of tunnelling amplitudes 
will break down. We refer to a pilot study [|Il|] on S 3 x IR for readers interested in this 
issue, and for an explanation of the relevance of the geometry T 3 x IR which is studied in 
this paper. This geometry allows us to find the instantons using the lattice approximation. 
For simplicity we restrict ourselves to SU (2) pure gauge theories. 



2 On the existence of continuum solutions 

The geometry T 3 x IR, in particular in a lattice formulation, can be seen as a limiting 
case of an asymmetric four torus [0,L] 3 x [0,T]. The only known solutions have constant 
curvature |12|] and hence cannot correspond to vacuum to vacuum tunnelling, furthermore 
their topological charge is at least 2. Actually, it can be proven rigorously [0, that for T 



finite, no regular charge 1 self-dual solutions can exist on a four-torus (we will illustrate 



this with our numerical results). As soon as we allow for twisted boundary conditions |14 
existence of minimal non-trivial topological charge instanton solutions can be proven. One 
distinguishes two cases, depending on the properties of the twist tensor n^ v G TL 2 . 

When ^e^xcrn^nxrj = lmod2, the topological charge is half-integer. The minimal 
action allowed by the topological bound is therefore 4n 2 , corresponding to topological 
charge 1/2. As twist is also well defined on the lattice [[HJ, and in the above situation (called 
non-orthogonal twist) does not allow for zero-action configurations, these instantons cannot 
"fall through the lattice" . Indeed, the index theorem predicts in this case 4 parameters (8 
x topological charge), which have to correspond to the position parameters. The charge 
1/2 instanton hence has fixed size and cannot shrink due to lattice artefacts. Impressively 
accurate results [|nj were obtained for this case using the well known cooling method [T7J 
to find a solution of the (lattice) equations of motion, whose smoothness and scaling with 
the lattice volume leaves no room to doubt it provides an accurate approximation to the 
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continuum solution with action 4tt 2 . In the continuum, existence of smooth non-trivial (but 



not necessarily self-dual) solutions was proven by Sedlacek [18|, whereas theorem 3.2.1. of 



ref. |T9j states that the moduli space of self-dual solutions is isomorphic with a four-torus. 



When ^e^xan^nxa = 0mod2, also called an orthogonal twist, there are "twist eat- 
ing" |1| configurations, i.e. configurations that have zero action and are compatible with 
twisted boundary conditions (see also |21|]). For SU(2), it is not too difficult to show that 
as long as n^ u ^ mod 2 for some /i and v, this twist eating configuration is unique , up 
to a global gauge transformation if a twist is introduced as in ref. ^] and multiplica- 
tion with elements of the center of the gauge group. With twisted boundary conditions as 



originally defined by 't Hooft [[14|, such a global gauge transformation would even change 
the boundary conditions, and as 5*0(3) bundles the twist eating configuration is unique 
(For SU(N) it can be proven [^| that out of the iV 4 center elements that can multiply the 
twist, only N 2 give rise to gauge inequivalent configurations). Under this condition it can 
be shown [jn| that there are instanton solutions with 8 parameters (its moduli space, when 
dividing out the trivial translation parameters, is even related to a K3 surface |[9L E0|) us- 



ing Taubes' |24| technique of glueing a localized instanton (with scale, position and global 
gauge parameters) to the "twist eating" flat connection (i.e. zero action configuration). As 
the latter is not invariant under global gauge transformations, the global gauge parameters 
of the localized instantons are genuine parameters of the moduli space (see also P5|]). 

The reason twisted boundary conditions are useful, is that at finite T there are no exact 
instantons on T 4 with periodic boundary conditions, but there are exact solutions for any 
non-trivial twist in the time direction. As T — ► oo these solutions are also solutions on 
T 3 x 1R. This comes about as follows. Since at T — > oo the action can only stay finite if 
for \t\ — > oo the energy density goes to zero, we deduce from a vanishing magnetic energy 
that up to a gauge 

Mx,t^ ±oo) = Cfa 3 /2L , (1) 

where Cf G [0,47r] (A = 0) parametrizes the vacuum or toron valley p6| , whose gauge 
invariant observables are best described by the Polyakov line expectation values 

Pi = iTr (p exp(* £ A{x, t)dxA = cos (d/2) , (2) 



(for the proper definiton in the presence of twist, see ref.||16||.) In the vacuum valley P, 



is space independent and the vanishing of the electric energy at t — ► ±oo also requires 
Pi (or Cf) to be asymptotically time independent. Instanton solutions on T 3 x R 
are hence characterized by the boundary conditions Cf at t — > ±oo. It is these general 
instantons that are physically relevant. It is not clear if solutions exist with arbitrary 
boundary values. Approaching T — > oo, by using periodic boundary conditions (which 
would impose Cf = Cf mod 47r up to a periodic gauge transformation) does not allow 
us to prove existence. As long as T is finite there are no solutions |L3| and the proof of 
non-existence breaks down as T — ► oo. On the other hand, with twist in the time direction, 
riQi = 1, even at T finite there is in the continuum an 8 parameter set of exact instanton 
solutions, which at T — > oo will correspond to Cf = (2n — Cf) mod47r (again up to a 
periodic gauge transformation). For localized instantons, asymptotically the field has to 
coincide with the unique flat connection, which fixes the possible values of Cf to 7r, but at 
the other extreme, as the instanton in the spatial direction extends up to the "boundary" 
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of the torus, the regions t —> +00 and t — > —00 no longer are connected, which will relax 
the fact that Pj = (Cf = ix). Although we have no proof, it is reasonable to assume that 
the 8 parameters for the instantons close to the maximal size are described by p, the 4 
position parameters and the 3 vacuum valley parameters Cf or C^). Note that for Pj — > 
as t — > ±00, the solution is both compatible with twisted and periodic boundary conditions 
at infinite T . In any case we have now learned that on T 3 x R (i.e. with free boundary 
conditions at t — > ±00) there are, at least 8 and at most 11 continuous parameters that 
describe the instanton solutions for vacuum to vacuum tunnelling. 



3 The lattice actions and cooling 

Let us start with discussing the standard Wilson action || 

S= E Tr ( 1 -0)= Y, W - Uv(x)U v (x + ftUt(x + 0)Ut(x)) , (3) 

where U^(x) are SU(2) group elements on the link that runs from x to x + //, the latter 
being the unit vector in the /1 direction. To derive the equations of motion, we observe 
that S depends on U^x) through the expression: 

S(U,(x))=Tt(1-U,(x)UI(x)) + Tt(1-UI(x)U,(x)) , (4) 

where 

(5) 

which is independent of U^x). Hence, S(e x U tl (x)) — S(U f j i (x)) = 0(X 2 ) for any Lie algebra 
element X, implies 

TtWUrWUliz) - U,(x)Ul(x))] = , (6) 

where Oi are the Pauli matrices. This is easily seen to imply that U^(x)W(x) is a multiple 
of the identity, and as is the sum of SU(2) matrices, it can be written as = a + ia- a, 
with a M G 1R 4 . If we define \\U^\\ = ^/a|, eq.© is seen to imply 

U,(x) = ±%(x)/\\%(x)\\ . (7) 

As we are only interested in stable solutions (i.e. local minima of the action), the plus sign 
in eq.(|7D is the relevant one. The process of iteratively finding the solution to the equations 
of motion is called cooling ||17|| , as in all cases it is devised such that the action is lowered 



after each iteration. The easiest is to simply choose U'Jjc) = U^{x) /\\U^{x)\\ since the fixed 
point of this iteration is clearly a solution to the equations of motion. An optimal way to 
sweep through the lattice is to divide for each \i the links U^(x) in two mutually exclusive 
checkerboard patterns IT^ such that all links on a particular pattern LT^ (i.e. for fixed % and 
/i) can be changed simultaneously, which is a well known trick to vectorize this procedure. 
At the cost of roughly a factor two in memory-use, vectorization is also achieved for the 
modified action we have considered so far for our numerical simulations: 



u x,fi,i/ * H> w x,[i,v \ 



(8) 
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The meaning of the parameter e will become clear in the next section. For ease of our 
numerical studies we have not considered modified single plaquette actions (see also the 
next section for a discussion on the adjoint and Manton actions). 



4 Lattice artefacts 

To calculate the effect of the discretization on the solutions of the equations of motion 
we first take a smooth continuum configuration (not necessarily a solution) A^x). For 
definiteness we put L = 1, and N s the number of lattice points in the spatial direction such 
that a = l/N s . We put this configuration on the lattice by defining: 

U^{x) = Pexp(i / A^x + sp,)ds) . (9) 

J 

The value of the plaquette thus corresponds to parallel transport around a square and 
can easily be proven to be given by |27| (D^ = + A^x) the covariant derivative in the 
fundamental representation) 

^ r (^/O = ^{e aD ^ e aDv ^ e~ aD,J -^ e~ aD '-'^) . (10) 
The proof simply amounts to observing that if A^x) = A^, i.e. A^ is space-time inde- 
pendent, then Tr = r Ii{e aA ^e aAv e~ aA ^e~ aAv ) and eq.flTUD is the only way to make 
this formula gauge invariant under arbitrary (i.e. x-dependent) gauge transformations. 
Using the Campbell-Baker-Hausdorff formula, eq.(^) can be expressed in terms of prod- 
ucts of covariant derivatives T>^ (in the adjoint representation) acting on the curvature 
Ffo, = [D^Du] = dA^ - dA u + [A„,A V ], e.g. V^F^ = [D„, D v \\. As the action 
involves a sum over all x, // and z/, things can be considerably simplified by computing, 
what we will call, the clover average 

Tr (jQ)) dover = j^Tr f J = ^Tr[e- aD »e~ aD " e aD »e aDv + e - aD » e aD » e aD » e - aDv 



= Tr[l + - ^ ((V,F^(x)Y + {V u F, v {x)f) + 

+^ ((^(^)) 2 + (PtF^)) 2 ) + \l%Fr{x)T%Fr{x) 
-hj)^V p F^(x)) 2 }} + C(a 10 ) + total derivative terms , (11) 

for which the multiple Campbell-Baker-Hausdorff expansion of eq.(fL0|) is required to 0(a 6 ), 
obtained with the aid of the symbolic manipulation program FORM |23. The clover 
average allows one to ignore many terms (all those odd in any of the indices) in evaluating 
the trace of the exponent. 

Equation ( |ll~D was also derived using the non-abelian Stokes formula (sq = 1) 

Ufatix) = Up{x)U v {x + pt)U^{x + u)Ul(x) = Pexp(a 2 J ds J dt T^ix + as fx + atv)) 

00 n [Si-! rl 

= 1 + ^2 II / / dU c^T^ix + asxjl + ativ) . . .a^T^ix + as n fi + at n v) , (12) 



n=l i=l 
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where T^iy) equals F^ v (y) up to the backtracking loop that connects y to x, or: 
V(s, t) = Pexp (a J A^(x + asfi)ds)Pexp(a J A^(x + asfx + atP)dt) , 

3~nv{% + as fx + atu) = V(s, t)F^ ll/ (x + asfx + atv)V*(s, t) . (13) 

To obtain the result of eq.(|ll|) one now expands T^vix + asfx + atu) around the point x, 
making use of the identity: 

d;d^ u (x)=v;v™F, u (x) . (14) 

Note that the ordering of the covariant derivatives in the r.h.s. of eq.([L4]) is essential. Also 
crucial is that the path ordering U(s,t) = Pexp( /* A{u)du) (where A{t) = e^A^x + te) for 
some unit vector e) is compatible with the covariant derivative, i.e. e^D^x + se)U(s,t) = 
= e^D^(x+te)W(s, t) (in this respect we have corrected the formula in ref. ||29|| ). Inserting 
the Taylor expansion of J r fJjU (x + asfi + atC') with respect to (s,t) in eq. fll2]) , gives the result 
of eq. ([□]). A very useful check is that the symmetry implied by U^ v (x) = (x), not 
explicit at intermediate steps of the calculation, is respected by the final result. 
Using eq.(|,jll]), one finds to 0(a 10 ) for the modified action S(e) 



+ 30 ( V l F ^ X W + (^ F ^( x )) 2 j + ^lF^x)VtF^x) - - A {V,V v F, v {x)) 2 }]. (15) 

Obviously, S(e = 1) corresponds to the Wilson action, and the sign of the leading lattice 
artefacts are simply reversed by changing the sign of e. Most of the numerical results were 
obtained for e — — 1, but e is useful in the initial cooling from a random configuration. By 
keeping e > as long as S > 8tt 2 , and only switching to e = —1 when S ~ 8tt 2 , we can 
avoid the solution to get stuck at higher topological charges. Once we set e = —1, we have 
yet to see an instanton fall through the lattice. We will come back to these issues when 
discussing the numerical results. Also note that, as Tr a d(t/) = |Tr(£7)| 2 — 1, one finds that 
the Wilson action in the adjoint representation (S^d) satisfies = 4S(e = 1) + 0(a 8 ) 
and does not allow us to change the sign of the a 6 term. The same holds for the Manton 
action |3(| which by definition agrees to 0(a 8 ) with the Wilson action. 

In the past, more complicated improved actions were considered || |31|, for which we 
will present the result similar to eq. flT5|) , as it allows us to predict whether or not they give 
rise to stable solutions ||. It also allows comparison with earlier results by Liischer and 
Weisz p]| obtained from a perturbative analysis. In the following, the coefficients in front 



of the Cj are to match with the definitions of ref. |]3T|. The averages (• • •) are similar to 
the clover average above, but include now also averaging over all orientations of the loops. 
After some algebra one finds: 



S({ Cl }) = £ Tr{c + 2d (i- ) + 4c 2 (i- J^) + ^c 3 (i- £1)} 



4 

--(c + 8d + 16c 2 + 8c 3 ) E Tr(F^(a;)) + a 6 (c 2 + |) £ Tr(V,F^(x)V u F uX (x)) 
+^(c + 20 Cl + 4c 2 - 4c 3 ) £ Tr(P M F^(x)) 2 + a 6 | £ Ti((V,F vX ) 2 ) + 0(a 8 ) . (16) 
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One can therefore achieve tree-level improvement by choosing |yj C0 + 8C1 + I6C2 + 8C3 = 1 



c + 20ci + 4c 2 — 4c 3 = and c 2 = c 3 = 0. Note that the condition c 2 + c 3 /3 = only 
applies off-shell, since on-shell J2x,uv\ Tr(P /J F AtA (x)T , !/ F jyA (x)) = 0. 

Iwasaki and Yoshie || considered cooling for the Symanzik improved action, that is 
Co = |, cx = — jz and c 2j3 = 0, for which the a 6 term vanishes. The a 8 term will have 
to be computed to settle stability. From equation (JTJp one sees that the a 6 term has a 
definite sign. This is no longer the case for the a 8 term. The same holds for the Symanzik 
improved action: 



a 4 „ , , a 8 r „a , x I 



x,/j,,v Z O 



. x 



- \{V,V u F, v {x)f + ^{VlF, v f}] + 0{a^) . (17) 

To decide in these cases if the lattice admits a stable solution (i.e. its action increases with 
decreasing p), one can compute the lattice action using explicitly the topological charge- 
one instanton solution with scale p. Eqs. ( |H||T6|) and ( |TTD are only valid as long as a C p 
because for p ~ a the expansion in powers of a no longer converges. For p C L to a good 
approximation we can substitute the infinite volume continuum instanton solution: 



AJx) = -i- 



IflV" 



with rj^ u the self-dual 't Hooft tensor |p3|| . When p ~ L the solution will of course be 
modified by the boundary effect. Substituting eq.(^) we find 

S(e) = ^{l--{a/pf-^^{a/pY + 0{a/pf} , 

Symanzik = 8tT 2 { 1 - ^ (a/ P ) 4 + 0{a/pf} , (19) 

we thus confirm the observation of Iwasaki and Yoshie || that the Symanzik tree-level 
improved action has no stable instanton solutions. Since S(e = 0) = 87r 2 {l + Jr(a/p) 4 + 
0(a/p) 6 } we predict even at e = the lattice to have stable solutions, which we have 
verified for the case with twisted boundary conditions in the time direction (see below). 

Iwasaki and Yoshie || also considered cooling for Wilson's choice [^] (W) of cq = 4.376, 
c\ = —0.252, c 2 = and c 3 = —0.17 and for (R) c = 9, c\ = —1 and c 2 = c 3 = 0. To 
0(a 8 ) these actions effectively correspond respectively to e = —2.704 and e = — 11, which 
for the case (W) we computed by substituting the continuum instanton solution. Indeed, 
they see stability up to 250 sweeps in both cases. 



5 Non-leading lattice artefact corrections 

In presenting eq.(^) we have replaced the sum over the lattice points by an integral and 
ignored the fact that on the lattice the equations of motion are modified. Both effects 



7 



turn out to be small, the first exponential in p/a, the other gives a correction to the 
expression for S(e) in eq.([T9|) proportional to e 2 (a/p) 4 (whereas the correction to Ssymanzik 
is proportional to (a/p) 8 , which also holds for S(e = 0)). 

We wish to compute J2 X f{ x ) = En f{ na )i f° r which we can use its Fourier decomposi- 
tion 

a 4 ]T/(x) = a 4 J2 J2e ik ™f(k) = a 4 N*N t £ /(^) = E / ^ xla f(x)d*x . 

(20) 

The terms with p ^ give the error one makes, when replacing the lattice sum by an 
integral. For a <C p <C L and f(x) = — ^Tr(F 2 (x+xo)) one finds explicitly (using eq . (|18|) ) 

-y E ^(^O* + *<>)) = 8tt 2 [1 + E 27rV(p/a) 2 co S (27rp • x /a)K 2 (2n\p\p/a)} 

x ^ v pG^ 4 \{0} 

= 8n 2 [l-87r 2 (p/a) 3/2 e~ 2np/a (l + 0(a/p))] , (21) 

(with K2 the modified Bessel function ||34j| ). Here we have taken xq to coincide with a 
point on the dual lattice, 2xq = a for all p, as this minimizes the action. 

To estimate the shift in the equations of motion due to the lattice artefacts, we again 
consider /)>a, such that the action can, in a good approximation, be given by 

= E J *x{--T*{F%{x)) + — Tr((V,F, u (x)) 2 )} + 0(a 4 ) , (22) 

which implies the equations of motion: 

YP»Fv* = ^ 2H , = + l -[F^V^ u ] + ^D\V V F VV ) . (23) 

As eq.(|22|) breaks the scale invariance, there will in general not be solutions close to eq. (|18|) . 
Variation with respect to p no longer leaves the action invariant. Still, since this variation 
corresponds to a near zero-mode, it makes sense to expect quasi-stability under cooling. 
The action changes only slowly in the direction of this near zero-mode but is predominantly 
lowered in those directions that leave the curvature square integrable and are spanned by 
the non-zero modes of the quadratic fluctuation operator for the action, which in the 
background gauge corresponds to 

M Xa = 5 Xa Vl + 2 a dF Xa . (24) 

If V is the projection operator on the normalizable non-zero modes of M. one has (at p = 1) 

A, = Af+ea 2 VM-lVH v ^Af+ea 2 A^ , 

W = (T^^' 1 W V) , (25) 

where H^(x) is evaluated by substituting for A^(x) the continuum solution A^(x) given 
in eq . (|l8l) . For convenience we introduced the quantities tp^'{x): 

.V2'2//„\ ' \ „0 „, _ - /A 2 ' 2 / / 



tpT 2> (x) = iJ2v a ^x»°a , ^T 2 \x) = iY,V a ^uCT a (x 2 -6x1) ' 

u,a u,a 

= tJ2v^x u a a (3x 2 -3xl-5xl) , (26) 



i/.a 
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which are eigenfunctions of the angular momentum operators L\ and J 2 (as defined in 
L\ = —^rj^Xpdv, J a = Lf + ad(^f)). To compute M~lVH u one can use for M~l the 
explicit expression |35| M,~l = VuX^xi^a^Via^^ ( m the gauge V^(VH^) = 0). The 



result, which can be verified by applying M.n U} is found to be 



M -i PH _ f lQ g(i + ^ 2 ) i i u (U) (x) 

M ^ FH » - ( 5 (1+* 2 ) 2 3(1 +X 2)3+ 10 (1 + X 2f " " 

21og(l + x 2 ) 6 + 3x 2 - x 4 (|, |) 
l 5x 8 (l + x 2 ) 2 15x 6 (l + x 2 ) 3 ^ M 1 
2(3 + 5a; 2 )log(l + a: 2 ) 6 + 13x 2 + 4x 4 . ,(§,§) 



+ ( 5x 8 (l+x 2 ) 2 5x6(1+ a; 2 ) 3 ^ ^ (x) ' (2?) 



The algebraic manipulation program Mathematica |36[ was useful in obtaining and checking 
these results. Despite its appearance, this result is regular at x —>■ 0. However, it contains 

a non-normalizable deformation (since ip^' 2 \x) /(l + x 2 ) = — A^(x)), which would make 
the action diverge and should be removed by projecting on normalizable deformations: 

4 1} = VM~JVH V = M-JVH V - { ^ ll tU r\x\ M^ = H, 



» ^ u v v 60(1 + a; 2 ) 2 rM K h ^ v ~ 11 5(1 + x 2 ) 3 ' 

(28) 

One easily verifies that Affl and Ai^A^ are both square integrable and orthogonal to the 



zero-mode dAf/dp^ = 2^1'^ (x) / (1 + x 2 ) 2 . 

We can now substitute A^ = Aff' + ea 2 A^ in eq.(^) (for later convenience evaluated 
for a different e) to obtain the shift in the action 

§(e) = S{e)-2eea 4 J d 4 xTr(A^{x)H^(x))+e 2 a 4 J d 4 xTr(A^(x)M fll/ A^(x)) + 0(a 6 ) , 

(29) 

where for the term linear in A^' we used the equations of motion for A^\ We also rein- 
troduced the p dependence using trivial dimensional arguments. Evaluating the integrals 
gives 

S(e) = 8n 2 {l - ! A 2 - [ 15e ~ - 12 + ggf! - lZgg! ]( ^)4 + £ f] . (go) 
KJ 1 5> ; L 210 2625 5250 J V V 

At e = e = — 1 the shift due to the modified equations of motion in the (a/ p) 4 term is 58%. 

Strictly speaking our expression for the p-dependence of the lattice action is only valid 
for p/a ~ 2 and p < L, since we are using the continuum infinite volume solution as the 
zero-order approximation. But even for p ~ L/2 it is not unreasonable to expect the order 
of magnitude of the corrections to be given by eqs.(^) and (|3"0|) . 



6 Numerical results and discussion 

This section discusses the numerical results obtained as described in section 3, mainly to 
illustrate the viability of our ideas. A more detailed and careful analysis will be left for a 
future publication. So far we have worked mainly on lattices of size N 3 x A^, with N s = 7 
or 8 and with N t = 3N S to 4N S . At e — —1 (see eq.(|])) we settle to an action near 87r 2 , and 
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we have seen stability for up to 6000 sweeps. The same is true for e = in the presence of 
a twist (but not without twist, where our configuration ultimately decays to the vacuum 
at £ — 0, but note that in that case there are no regular instanton solutions). Apart from 
the total action we compute separately the sum over the n x m plaquettes, denoted by 
S nxm and averaging over the two orientations if n ^ m. S nxm is normalized by dividing 
by 8n 2 n 2 m 2 , such that for an infinite lattice and p/a — > oo, S nxm — > 1. When we perform 
cooling with the action of eq.® we should take = A^> + ea 2 A^ in calculating S nxm . 
Eq.(|TT]) for Sixi easily leads to the general result for S nxm by inserting for each index p(y) 
a factor n(m). Together with eq.(p9|) one deduces, that to C(a 6 ) 



1 



(n 2 + m 2 ) 



-aa 



m 2 n 2 (3i 



4 i 4 

m + n 



(3 2 + 



n 2 + m 2 



e 7 - e 2 S a 4 , (31) 



up to the discretization error implied by eq. (|21|) , which for the lattices we are considering 
can be estimated to be not bigger than 10~ 6 . This formula holds for sufficiently smooth 
configurations, i.e. aa 2 <C 1, even if the configuration has non- vanishing action over the 
entire spatial volume. It is these configurations that are of interest to us and which deviate 
considerably from localized instantons (eq.(|T8|)) for which p <C L. From eqs. ( |T9| , |30| ) we 
easily deduce for those localized instantons the results: 



aa 



U a -: 

5 p 



Pia 4 



29 a N 
630>' 



63 V 



7a 



284 a, 
2625 V 



5a 4 



179 
5250 V 



(32) 

From the numerical results we have obtained Si x i, £1x2, 5*2x2 and Si X 3 on two lattices of 
size respectively 7 3 x 21 and 8 3 x 24, for e = and e = —1 with a twist n 0i = (1, 1, 1) (see 
table I). From these we extract the coefficients in eq. (|3lD , whose values are summarized in 
table II (the error due to neglecting the 0(a 6 ) term is of the order of (n 6 + m 6 )(ao 2 ) 3 ). 



N s 




e 




5*1x2 


5*2x2 


5lx3 


7 3 


x 21 


-1 


0.982591 


0.957050 


0.928823 


0.918105 


7 3 


x 21 





0.982287 


0.956437 


0.927908 


0.917109 


8 3 


x 24 


-1 


0.986720 


0.967122 


0.945887 


0.936619 


8 3 


x 24 





0.986529 


0.966736 


0.945310 


0.935976 



Table I: Numerical results obtained by cooling with S(e) and twist noi = (1, 1, 1). 



N s xN t 


aa 2 


Not 2 


No? 


7/a 2 


S/a 2 


7 3 x 21 


0.01761 


0.96 


0.63 


0.66 


0.32 


8 3 x 24 


0.01340 


1.01 


0.64 


0.71 


0.35 


p « L 


0.2(a/p) 2 


1.15 


0.79 


2.70 


0.85 



Table II: Coefficients appearing in eq.(|31~D extracted from the numerical 
results in table I, using 5i x i, 52 X 2 and 5i X 2 (the latter at e — — 1 only). 
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It is interesting to analyse the untwisted case in more detail to illustrate the difficulty in 
having self-dual solutions at finite T. In figure 1 we plot (a) the total electric and magnetic 
energies Se,b{^)i (b) the Polyakov-line P\(t) through two particular points x and similarly 
for P2,z{t) m ( c ,d). We see two features that are intimately related. First, where 8 Bit) — > 0, 
the electric energy £e(£) — ► const. Second for the same t values where this occurs Cj(i) 
(Pi(t) = cos(Cj(t)/2)) is linear in t and x independent. These are precisely the equations 
of motion when restricting to the vacuum valley. Classically motion on this valley, which 
itself has the geometry of a three-torus, is free. On the lattice this motion is described by 
the action 

£4A?(l-cos(^±±L^)) . (33) 
One easily checks that the values of C&(t + 1) — Ck(t) obtained from figs.l(b,c,d) quite 



accurately reproduce through eq.(|33f) the value for £e(t). Clearly the electric tail destroys 
the self-duality. Suppose that at T — > 00 the solution describes tunnelling from Cf to Cf 
and Cf 7^ Cf , then at finite T the periodic boundary conditions force Ci(t) to linearly 
interpolate between Cf and Cf over a time T — T , if T is the time interval over which 
8b(£) 7^ 0. Thus the action, even in the continuum, would be bigger than 871" 2 by a number 
proportional to 1/(T — T ) except when there are solutions with Cf = Cf for T — > 00. 
These can certainly not be excluded, in particular as Cf = Cf = it is compatible with a 
twist riQi = 1, but if these are very localized instantons, the lattice artefacts might make 
its action so big, that the lattice will prefer the least localized solutions with Cf 7^ Cf. If 
T is not big enough the lattice will find a compromise between these two cases. There are 
indications that the largest instanton prefers Cf 7^ Cf and from the numerical results with 
twist n 0i = (1, 1, 1) presented in fig. 2, the preferred values seem to be such that two of the 
Cj go from to 2tt and one goes from 2tt/3 to 47r/3. We compare (after appropriate scaling 
with N s ) for N s = 7 and 8, using over- improved cooling at e — — 1, in fig. 2a the electric 
and magnetic energy profiles, and in figs.2b-d the values of Cj = 2acos(Pj), at the spatial 
lattice point with maximal energy (to be precise, with maximal Ef). From this we deduce 
8e = 8b to a high accuracy, consistent with self-duality, and the excellent scaling with N s . 
The results in figure 2 are obtained after roughly 6600 cooling sweeps, which is necessary 
since the dependence of the lattice action on Cf is rather weak (at e = too weak to 
observe) and the configuration only slowly reaches the minimum of the lattice action. We 
have verified that the approach to this minimum is exponential, as is illustrated in fig. 3, 
where we plot the total action and the maximum of 8s(t) (i.e. £#(0)) as a function of the 
number of cooling sweeps. We see indeed that the maximal energy along the tunnelling 
path decreases under cooling, which is mainly due to the increasing size, as otherwise the 
action should depend more strongly on the number of cooling sweeps. (For the Wilson 
action one sees a dramatic increase of £#(0) under cooling, until the action suddenly drops 
to zero.) With boundary conditions that fix the link variables at t — and t — T to the 
vacuum configurations, the approach to the minimum action is much faster. 

Elsewhere we will publish a more detailed analysis of the scaling properties, as well as 
testing our conjecture to be able to find a sphaleron. Also numerical results with fixed 
boundary conditions, that allow us to investigate if solutions exist for arbitrary Cf will 
be presented elsewhere. This paper mainly served the purpose to describe the formalism, 
and demonstrate the large amount of control obtained in this way in studying instanton 
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solutions on a torus. 

It would be interesting to repeat this analysis for the 2-dimensional 0(3) model, for 
which the instantons on a torus are exactly known |j37| , in the light of the "perfect" lattice 



action recently considered by Hasenfratz and Niedermayer ||38|| . But as we have shown, 
appropriate deviations from a "perfect" action can be quite helpful. 

Finally, over-improvement might be an efficient tool to measure the topological suscep- 
tibility, as the action to generate a statistical ensemble need not be the same as the one 
used to measure the topological charge. 
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Figure captions 



Figure 1: Numerical results (after scaling appropriately with N s ) for the case 
of an 8 3 x 24 lattice without twist, obtained from over- improved cooling at 
e = — 1. In (a) the electric triangles) and magnetic (£b(£) squares) 

energies are plotted. In the upper part of this figure the tails are plotted 
at an enlarged scale. In (b-d) are plotted Ci(t) = 2acos(Pi(^)) through two 
distinct spatial points on the lattice. 

Figure 2: Numerical results (after scaling appropriately with N s ) for the 
cases of a 7 3 x 21 (squares) and an 8 3 x 24 (triangles) lattice with twist 
rioi = (1, 1, 1), obtained from over- improved cooling at e = —1. Figure (a) 
contains four data sets. Two for with the above mentioned symbols 

and two (crosses for N s — 7 and stars for N s = 8) for Figures (b-d) 

exhibit C;(t), through the spatial lattice point with maximal E\ at t — 0. 

Figure 3: The history of the action S(e = — 1) and the maximal magnetic 
energy £ B (t = 0) as a function of the number of cooling sweeps for an 8 3 x 24 
lattice with twist n 0i = (1,1,1), together with their exponential fits. The 
short lines on the right indicate the asymptotic values following from these 
fits. 
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